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Abstract 

We present the MBR (Minimum Bias Rockefeller) Monte Carlo simulation of (anti) 
proton-proton interactions and its implementation in the PYTHIA8 event generator. We 
discuss the total, elastic, and total-inelastic cross sections, and three contributions from 
diffraction dissociation processes that contribute to the latter: single diffraction, double 
diffraction, and central diffraction or double-Pomeron exchange. The event generation 
follows a renormalized-Regge-theory model, successfully tested using CDF data. Based 
on the MBR-enhanced pythiaS simulation, we present cross-section predictions for the 
LHC and beyond, up to collision energies of 50 TeV. 



1 Introduction 

The MBR (Minimum Bias Rockefeller) Monte Carlo (MC) simulation is an event generator 
addressing the contributions of three diffraction-dissociation processes to the total-inelastic 
pp cross sectioiiQ single-diffraction dissociation or single dissociation (SD), in which one 
of the incoming protons dissociates, double-diffraction dissociation or double dissociation 
(DD), in which both protons dissociate, and central dissociation (CD) or double-Pomeron 
exchange (DPE), where neither proton dissociates. These processes are tabulated below, 

SD pp — > Xp 

or pp — )■ pY 

DD pp XY 

CD (DPE) pp pXp, 

where X and Y represent diffractively dissociated protons. Schematic diagrams are shown in 
Fig. [T| along with diagrams for the total and elastic cross sections, which are also simulated 
in MBR. 

MBR predicts the energy dependence of the total, elastic, and total-inelastic pp cross 
sections, and fully simulates the above three diffractive components of the total-inelastic 
cross section. The diffractive-event generation is based on a phenomenological renormalized- 
Regge-theory model [Tj, originally developed for the CDF experiment. We have implemented 
MBR in PYTHIAS^, where it can be activated with the flag: Dif fraction :PomFlux = 5. 

^Although the original mbr addresses several hadron-hadron collisions, including pp, we will be assuming 
pp collisions throughout this paper for simplicity, except as explicitly stated. 
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Figure 1: Schematic diagrams of soft pp processes addressed in mbr: (a) total cross sec- 
tion, (b) elastic scattering, (c)-(d) single diffraction or single dissociation (SD), (e) double 
diffraction or double dissociation (DD), and (f) central diffraction (CD) or double-Pomeron 
exchange (DPE). 

This paper is organized as follows: in Sec. [2] we present the formulae used to calculate the 
total and elastic cross sections, as well as the diffractive and non-diffractive components of 
the total-inelastic cross section; in Sec. |3]we describe details of the generation of diffractive 
events; and in Sec. |4]we summarize the pythiaS "steering cards" for employing mbr in the 
simulation. 



2 Cross sections 

The total cross section (o"tot) is the sum of the total-elastic (dei) and total-inelastic (o"inei) 



cross sections. The cxtot and are calculated in Sec. |2.1[ The total-inelastic cross section 
(o"inei) is obtained as cxinei = ctot — cTei- The cxinei receives contributions from the diffractive 
components (SD, DD, CD or DPE) and from the non-diffractive (ND) cross section, defined 
as: 

0"ND = (cTtot - 0-cl) - (2crsD + (^DD + CTCd), (l) 

where is the cross section for either pp — > Xp or pp — pY, assumed to be equal. The 
diffractive cross sections are calculated in Sec. 12.21 

2.1 Total, elastic, and total-inelastic cross sections 

The cr^^t^{s) cross sections at a pp center-of-mass-energy ^/s are calculated as follows: 
• i6.79s0.i04 ^ 60.81s"°-32 ^ 31.68s-°-^^ for < 1.8 TeV, 

(2) 

for > 1.8 TeV, ^ ^ 



p^p 

'^tot — \ ^CDF 



"tot ^ so 



2 ^ CDF ^ 2 



In ^ - In 



The term for y/s < 1.8 TeV, where the (±) denotes (p), is obtained from a global 
Regge-theory fit to pre-LHC data on p'^p, K^p and n^p cross sections [3], while that for 
> 1.8 TeV is a prediction of a model based on a saturated Froissart bound[4]. The 
latter, which we normalize to the CDF measurement of dtot at V s^°^ = 1.8 TeV, cr^^^ = 
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80.03 ± 2.24 mb, depends on two parameters: the energy at which the saturation occurs, 
■y^sp = '^fl GeV, and the energy-scale parameter, sq, for which we use sq = (3.7 ± 1.5) GeV^, 
divided by {hc)^ ~ 0.389 GeV^mb to obtain the cross section in Eq. (^in mb. 

The elastic cross section, cTgi ^, is calculated using atot from Eq. multiplied by the 
elastic-to-total cross-section ratio, cTd/cTtot, obtained from the global Regge fit of Ref. [3]. 
The total inelastic cross section is calculated as o"inci = crtot — <^e\- 

The energy dependences of cJtot, Cei and cxinei are shown in Fig. [2l and cross-section values 
at a/s=0.3, 0.9, 1.96, 2.76, 7, 8 and 14 TeV are presented in Tab. [ij 



2.2 Diffract ive cross sections 



Cross sections for SD, DD and CD (or DPE) are calculated using a phenomenological model 
discussed in detail in Ref. p]. Differential cross sections are expressed in terms of the 
Pomeron (iP) trajectory, a(t) = 1 + e + oit = 1.104 + 0.25 (GeV^^) • t, the Pomeron-proton 
coupling, f3{t), and the ratio of the triple- J' to the iP-proton couplings, n = g(t)/f3{0). For 
sufficiently large rapidity gaps {Ay > 3), for which iP-exchange dominates, the cross sections 
may be written as. 



dtdAy 
dtdAydyo 

d'^O'DPE 

dtidt2dAydyc 



gapV 
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,2[a(i)-l]As/ 
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167r 
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where t is the square of the four-momentum-transfer at the proton vertex and Ay is the 
rapidity gap width. The variable ?/o iii Eq. (|4]) is the center of the rapidity gap. In Eq. ([s]), 
the subscript i = 1,2 enumerates Pomerons in the DPE event. Ay = Ayi + Ay 2 is the 
total (sum of two gaps) rapidity-gap width in the event, and yc is the center in r] of the 
centrally-produced hadronic system. 

Eqs. ([3]) and ^ are equivalent to those of standard- Regge theory, as ^, the fractional 
forward-momentum-loss of the surviving proton (forward momentum carried by IP), is re- 
lated to the rapidity gap by ^ = 
^DD = MlM^/{s ■ So), where 
(DD) events. For DD events, yo = 

The Pomeron-proton coupling 



where /3(0) = 4.0728 Vmb 
F\t) = 



e The variable ^ is defined as ^sd = M'^/s and 
{Mf, M|) are the masses of dissociated systems in SD 
i ln(M|/M2), and for DPE ^ = = M^s. 
, (3{t), is given by: 

/3\t) = (3'{0)F\t), 
-1 



(6) 



6.566 GeV ^ and F{t) is the proton form factor from Ref. [5]: 

2 

(7) 



Ami ~ 2 



it 



Aml-t 



1 



1 



t 

0.71 



The right-hand side of Eq. ([7| is a double-exponential approximation of F'^it), with oi =0.9, 
02=0.1, hi =4.6 GeV"^, and 62 =0.6 GeV~^. The term in curly brackets in Eqs.([3])-(|5]) is 
the ]P-p total cross section at the reduced IP-p collision energy squared, s' = s ■ e ^ . The 
parameter k is set to k = 0.17 [9J, and /t/3^(0) = ctq, where Gq defines the total Pomeron- 



proton cross section at an energy-squared value of Sq 
7.249 GeV~^ 



1 GeV , is set to ctq = 2.82 mb or 
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Table 1: Cross sections in mb of the processes contributing to the Minimum-Bias sample, 



calculated as explained in Sec. 2.1 and Sec. 2.2, for selected values of ^/s. 



Pointing again to Eqs.([3|-([5| the renormalization parameter, A^gap(s), is defined as Ngg,p(s) 
mm(l, /), where / is the integral of the term in square brackets, corresponding to a Pomeron 
flux. The integral is calculated over all phase space in tj and in rjo (DD) or rjc (DPE) for 
Ay > 2.3. The renormalization procedure we use, which is based on interpreting the Pomeron 
flux as a (diffractive) gap-formation probability, yields predictions which are in very good 
agreement with diffractive measurements at CDF [6], [3, [S]. 

The cross-section formulae are used to generate events with large (diffractive) rapidity 
gaps. Ay. We suppress cross sections at small value of Ay by multiplying Eqs. (|3])-([5]) by: 



^"2 



1 + erf 



Ay - Ays \ 



(8) 



where erf is the error function [10], centered by default at Ays = 2 with a width of as = 0.5. 
For SD, this choice corresponds to suppressing events above the coherence limit < 0.135) 
[TT]. For DD, the choice of Ays = 2 is somewhat arbitrary, because of the difficulty of 
unambiguously distinguishing a low-A|/ DD event from a ND event with a rapidity gap from 
(exponentially suppressed) fluctuations. Changing Ays introduces event migrations between 
DD and ND samples. Such migrations have no impact on the total-inelastic cross section, 
as can be inferred from Eq. ([T]). For DPE, the suppression is applied to the total-gap width. 
Ay = Ay, + Ay 2. 

Fig. [3] presents the energy dependence of diffractive cross sections calculated as explained 
above. The non-diffractive cross section, given by Eq. Q is also shown. Cross-section values 
for v^=0.3, 0.9, 1.96, 2.76, 7, 8 and 14 TeV are tabulated in Tab. Q 



3 Event generation 

The four-momenta of particles produced in the event are generated as described below, in 



Sees. 3.1 3.3 The results of the simulation for SD and DD at ^/s = 7 TeV are presented 
in Fig. |4| which shows differential cross sections as a function of rapidity-gap width. Ay, com- 
pared to predictions of PYTHIA8-4c (rescaled Schiiler & Sjostrand model [2j, Dif fraction :PomFlux 
= 1). The distribution of da/dAy vs. the total rapidity-gap- width. Ay = Ay, + Ay2, and 
vs. the width of an individual gap, Ayi, for CD at ^/s = 7 TeV is shown in Fig. ^ 
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Figure 4: Differential cross sections as a function of Ay for (a) SD and (b) DD at ^/s = 
7 TeV, compared to predictions of pythia8-4C (rescaled Scliiiler & Sjostrand model[2], 
Dif fraction :PomFlux=l). 
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Figure 5: Differential cross sections for CD (DPE) at ^/s = 7 TeV as a function of: (a) 
total-gap width {Ay = Ayi + Ay2), and (b) single-gap width {Ayi). 



3.1 Single-difFr active events 

Events are generated by first choosing the rapidity-gap width, Ay, according to Eq. (|3| 
integrated over t: 



da 



SD ,Ay 



dAy 



ai 



+ 



hi + 2a' Ay 62 + 2a' Ay 



The range of the generation is defined by Ay„ 



Mq =MBRm2Min. The term: 



S 



^-S. (9) 

and Ay max = — InMg/s, where 



1 r ( Ay - MBRdyniinSD\ 
- 1 + er / 

2 [ V MBRdyminSigSD / 



(10) 
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is added to suppress events at low values of Ay, as explained in Sec. 2^ 
A value of t is then chosen according to: 



da 



SD 



dt 



F\t) 



where F'^{t) is given by Eq. (7) and the integration is performed up to 

-Ay 



with 



^ = e~^^. The diffractive mass is calculated as M = The four-momenta of the outgoing 
proton and the dissociated mass system are calculated using Mandelstam variables for a 
two-body scattering process, as implemented in PYTHIAS for other Dif fraction :PomFlux 
options. 



3.2 Double-difFractive events 

Events are generated by first choosing the rapidity-gap width according to Eq. Q integrated 
over t. Eq. (|4]) is divergent as Ay — )■ 0. In order to remove the divergence, the integration 
over t is performed within the limits from tmin = — e^^ to tmax = —e~^^. Then, Ay is chosen 
from the distribution: 

dAy 2a' Ay \ / ' ^ ' 

and the range of the generation is defined by Aymin = and Aymax = —^^Mq/{sso), where 
Mq =MBRm2Min and Sq = 1 GeV^. To further suppress events at low values of Ay the term: 



\ MBRdyminSigDD /J 



is used as explained in Sec. 2.2 



The variable t is chosen according to: 

^ ~ e'-'^y\ (14) 

in the range from tmin = — e^^ to tmax = — e~^^. 

Then, the center of the rapidity gap, yo, is selected uniformly within the limits: 

-K'"i-^^)<-<K'"i-^^)' 

and the diffractive masses are calculated as: 

= Vs ■ e-^y-yo, (16) 

M| = Vs ■ e-^y+yo. (17) 

The four-momenta of the outgoing dissociated mass systems are calculated using Man- 
delstam variables for a two-body scattering process, as implemented in pythiaS for other 
options of Dif fraction iPomFlux . 
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3.3 Central-difFractive (DPE) events 

Events are generated by first choosing the total rapidity gap width, Ay, according to Eq. ([s]), 
integrated over ti and 

— - e^^y / dyo f--U- S^S2, (18) 



dAy 
where: 



A2//2+yo 



f = I 1 ri9) 

•'^ \bi + a'Ay±2a'yo b2 + a'Ay ±2a'yoJ ' ^ ' 

and the integration is performed from Aymin = to Aymax = — In Mq / s, where Mq =MBRm2Min. 
For events at low values of Ay we suppress individual gaps with the factor: 

5=i[l + er/(^^-™'"^'^"^""°(_^)1. (20) 
2 [ VMBRdyniinSigCD/V2/J 

Then, the direction of the centrally-produced hadronic system, y^, is selected uniformly 
within the region: 

- 2 ^'^y ~ ^^™") <y<^< 2 ~ ' ^^^^ 

and rapidity gaps corresponding to each of the two Pomerons are calculated as: 

Ayi = Ay 12 + 2/0, (22) 

Ay2 = Ay/2 - y^. (23) 

The four-momentum transfers squared at each proton vertex, ti and are generated ac- 
cording to: 

'^^^F\u)-e^'^'^y^'\ (24) 

up to tmax,i = —m-l ■ whcrc = e~^^* and i = 1,2. Then, the pt and Pz of outgoing 

protons are calculated as ^ = (1 — ^j)|tj| — m^^f and \pz,i\ = p{^ — ^i), where p = ^ s/A — 
is the incoming proton momentum. 

Finally, the four-momentum of the hadronic system is calculated from the sum of the 
four-momenta of the Pomerons, each calculated as a difference between the incoming and 
outgoing proton four- vectors. 



4 MBR implementation in PYTHIA8 

The MBR generation is activated with the following flag: Dif fraction iPomFlux = 5. The 
simulation works for the pp, pp and pp scattering, and the beam setup is checked in the 
PYTHIA8 initialization step. It is assumed that the user will veto the event generation, after 
Pythia: : initO = false is returned for other beam particles. 

The simulation of the CD (DPE) process is implemented in pythiaS for the flrst time. It 
is activated with the flag SoftQCD : centralDiffractive = on. The corresponding-process 
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number is set to 106, naturally extending the list of soft QCD processed In the event record, 
the outgoing protons' information is written in rows 3 and 4, and the centrally-dissociated 
hadronic system occupies row 5 (represented by the rho_diff pseudo-particle). 

Below, we present the default values of parameters, used when the mbr simulation is 
activated: 



• the parameters of the Pomeron trajectory, a{t) = 1 + e + a't: 



Diffraction :MBRepsilon = 0.104 
Diffraction :MBRalpha =0.25 

the Pomeron-proton coupling, /3(0) (GeV~^), and the total Pomeron-proton cross sec- 



tion, ctq (mb), see Sec. 2.2 



Dif fraction :MBRbetaO = 6.566 
Diffraction :MBRsigmaO =2.82 

the lowest mass-squared value of the dissociated system, Mq, used to evaluate the 
highest allowed rapidity gap width, Ai/max, see Sec. [3] 



Diffraction :MBRm2Min = 1.5 



the minimum width of the rapidity gap used in the calculation of Ngap{s) (flux renor- 



malization. Sec. 2.2) 



Diffraction iMBRdyminSDf lux =2.3 
Diffraction iMBRdyminDDf lux =2.3 
Dif fraction iMBRdyminCDf lux =2.3 

the parameters Ays and as, used for the cross-section suppression at low Ay (non- 
diffractive region), see Eq. (Is]): 



Diffraction :MBRdyminSD =2.0 
Diffraction iMBRdyminDD =2.0 
Diffraction iMBRdyminCD =2.0 



Dif fraction iMBRdyminSigSD =0.5 
Diffraction iMBRdyminSigDD =0.5 
Diffraction iMBRdyminSigCD =0.5 

In addition, if the option Dif fraction :PomFlux = 5 is set, no dampening of the diffrac- 
tive cross sections is performed (no effect of SigmaDiff ractive : dampen = on); however, 
the user can still set his/her own cross sections by hand (SigmaTotal : setOwn= on). 

^103 and 104 for SD with proton dissociation along the positive or negative z-axis, respectively, and 105 
for DD. 



9 



5 Summary 

We present the mbr (Minimum Bias Rockefeller) Monte Carlo simulation, which has been 
tested using CDF data, and discuss its implementation in the pythiaS generator. The 
double-Pomeron-exchange process is included in pythiaS for the first time. The simulation 
is designed to apply to all Minimum Bias processes in the LHC energy range. 
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